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Numerical simulations of Blandford-Znajek energy extraction at high spin have revealed that held 
lines tend to bunch near the poles of the event horizon. We show that this behavior can be derived 
analytically from the assumption of fixed functional dependence of current and held line rotation on 
magnetic hux. The argument relies crucially on the existence of the Znajek condition, which offers 
non-trivial information about the helds on the horizon without requiring a full force-free solution. 
We also provide some new analytic expressions for the parabolic held configuration. 


I. INTRODUCTION 

A leading candidate mechanism to power relativistic 
jets from active galactic nuclei (AGN) is the Blandford- 
Znajek (BZ) process [T], in which energy is extracted from 
a spinning black hole via its plasma magnetosphere. In 
light of the large observed variation in jet properties, it is 
of interest to explore the detailed dependence of magne¬ 
tosphere structure on system parameters. One particular 
variation was noticed by [2] (TNM), who found that as 
the black hole spin is increased, the magnetic flux on the 
event horizon tends to concentrate near the pole, a fea¬ 
ture they described as bunching of field lines. In this note 
we will show how the bunching can be derived by a simple 
analytic argument from an assumption about the scaling 
of the current and field line rotation with spin. The ar¬ 
gument makes no reference to a particular choice of field 
line geometry, and we thereby show that the bunching is 
a general phenomenon, not limited to the particular ge¬ 
ometries considered by TNM. We illustrate the high-spin 
bunching for the three cases where approximate analyti¬ 
cal solutions are known, the radial mm, parabolic mm, 
and hyperbolic 0E3 field configurations. 


II. GENERAL MAGNETOSPHERE 
STRUCTURE 

Following BZ and TNM we assume that the plasma 
inertia is negligible, so that the magnetosphere is force- 
free 0 A stationary, axisymmetric force-free field in a 
spinning black hole background is fully characterized by 
the flux function if(r, 9), the polar current I(if), and the 
angular velocity of field lines 0(^)1^] These functions 
encode the energy flux by 

1=2/ dipl(if)fl(if), (1) 


where we assume a reflection-symmetric magnetosphere 
and define if on the northern hemisphere, with if zero 
at the pole and monotonically increasing until it reaches 
its largest value if* on the equator. (The factor of 2 ac¬ 
counts for the southern hemisphere.) Three families of 
energy-extracting approximate solutions are known (ra¬ 
dial, parabolic, and hyperbolic) and in all cases the cur¬ 
rent and angular velocity take the form 

i(if) = n H i(if), n = n H o(if), ( 2 ) 

where Qh is the horizon angular velocity and X and O 
do not depend on the spin a. This leads to the basic 
prediction of the BZ model that the power scales as the 
spin squared, 
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This result was derived analytically for small spin, but 
TNM has shown that for generic held configurations, it 
conti nues to hold for all but the highest (a > 0.99 M) 
spins T In fact TNM find more: not just the total power 
but also the detailed functional forms in Eqs. © carry 
over from the linearized theory to large spin. In par¬ 
ticular, their Figs. 5 and 6 demonstrate that I (if) and 
0(if) vary by no more than 10% over the entire range 
of spins 0 < a/M < 0.9999. We will take this result as 
our starting point, and for the remainder of the paper 
we will assume that Eqs. © hold exactly for any spin, 
with T(if) and 0(if) independent of the spin parameter 
a. We will also assume that the total magnetic flux if* is 
independent of spin. We now show how the bunching of 
held lines can be recovered straightforwardly from these 
assumptions. 
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1 Keviews ol force-tree electrodynamics, including the discussion 
of the quantities relevant to this work, may be found in [318]. 

2 We will use the conventions of TNM. 


3 As TNM show, it important to write the linearized result as , 
rather than (say) a * 1 2 , in order for the scaling to carry over to high 
spin. 

4 To define the notion of the “same” field configuration at different 
spins, TNM use the same initial data at each (r, 0) in Boyer- 
Lindquist coordinates. 
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A. Bunching of field lines 


A. Radial 


Even in possession of the detailed functional forms of 
I (il’) and (VO ) determining a complete force-free solu¬ 
tion is a daunting task, requiring the solution of a second- 
order, non-linear partial differential equation for ip. How¬ 
ever, for a Kerr black hole this equation has the remark¬ 
able property that, when evaluated on the horizon, only 
derivatives tangential to the horizon appear, and further¬ 
more the equation can be integrated once, so that it be¬ 
comes a first-order ordinary differential equation on the 
horizon. The result is the so-called Znajek condition m 
and e.g. 0), 
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where 1 Ih = a/ (tjj + a 2 ) is the horizon angular velocity, 
rjj = M + \JM 2 — a 2 is the horizon radius, and every¬ 
thing is evaluated on the horizon. We can integrate and 
then exponentiate Eq. © to find 

m = m, ( 5 ) 


where we defined 


/(VO = exp 
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/ 


and 


( 6 ) 


A(6) = e anHCOsB tanf. (7) 

Eq. ([6]) defines f(ip) up to an overall constant to be fixed, 
after solving Eq. ([5]), by demanding that the maximum 
value of ip be ip*. Formally, the solution to Eq. 0 is 

ip H (0) = f^im] = r 1 [e af2ff cose tan |] , (8) 


where the notation ip = ipu serves as a reminder that this 
formula holds on the horizon. Given the assumptions 
that X(ip), 0(ip), and ip* are independent of the spin, 
we know that /, and hence / -1 , is likewise independent 
of the spin. Thus the only dependence of the horizon 
flux on the spin is through the factor e aQ Hcos6 j n ^-j le 
argument of / _1 . This factor is monotonically decreasing 
from its maximum e a H at the pole to its minimum 1 at 
the equator, and will therefore always tend to increase 
the proportion of magnetic flux near the pole as the spin 
is increased. This is the bunching of field lines. 


III. SPECIFIC MAGNETIC GEOMETRIES 


We now consider force-free solutions with radial, 
parabolic, and hyperbolic geometries. In each case we 
use the current and angular velocity functions appro¬ 
priate to a normalization of ip(n/2) = ip* = 1. Fac¬ 
tors of ip* may be reinstated by scaling ip —> ip ftp*, 
I —» I/ip* and 17 —► 17. E.g., equation 0 becomes 
I (ip) = 2i:Q.(ip) ip (2 — ip/ip*). 


In their original paper BZ found an approximate solu¬ 
tion with radial field lines in a “split monopole” configu¬ 
ration. For ip* = 1 the current and angular velocity are 
given by 


I(ip) = 2nSl(ip) ip(2 - ip), (9) 

n (ip) = (io) 


which satisfy our assumptions with X = Trip(2 — ip) and 
O = 1/2. These are depicted in Fig. 1(a) which can 
be matched directly to Fig. 4 of TNM. The integral in 
Eq. 0 becomes 


/(VO = C. 


y> 

2 — ip’ 


( 11 ) 


for some constant C. We can then solve Eq. 0 on the 
horizon to learn that ip = 2A 2 / (C 2 + A 2 ). Requiring 
ip(w/2) = ip* = 1 then fixes C 2 = 1, and hence 


ip H (0) 


2 A 2 (9) 

1 + A 2 (6)' 


( 12 ) 


To compare directly to TNM we plot the radial mag¬ 
netic field B r = dgip/y/—g at the event horizon r = ru 
(Fig. 1(a)) where y/—g = (r 2 + a 2 cos 2 9) sin 9. Com¬ 
paring with their Fig. 7, we see excellent agreement for 
all curves a < 0.9, with TNM seeing more bunching for 
larger spins. This discrepancy can be explained by the 
fact that our assumption of spin-independent X(ip) and 
0(ip) disagrees more with TNM at higher spins. 


B. Parabolic 


BZ also found an approximate solution with parabolic 
field lines, whose current and angular velocity are 


I(lp) = 47T fl(ip) ip, 

, , __ (1 -VQ(m + 1)[m- (1-V»)ln2] 

H u 2 + (1 — ip) u — (1 — ip) 2 (u + 1) ln2’ 


(13) 

(14) 


where u = E[(l — ■0)ln4] and y = F(x) is the prod¬ 
uct logarithm, which is defined by the principal solu¬ 
tion of x = ye v ^\ Thus this configuration also has spin- 
independent X(ip) and 0(ip). 

Following the same steps as before, we find that 


m = c 


(1 — ip) In 2 


- 1 , 


(15) 


5 The angular velocity is normally given in terms of coordinates 
rather than as a function of ip. As far as we know, Eq. 01 is a 
new expression. 
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and for the horizon flux function that 


and 


ip H (d) 


A 1 2 3 4 5 6 (0) In 2 + ln[l + A 2 (6)] 
[1 + A 2 {6)] In 2 


(16) 




4A 2 (6>) 

[i+^wr 


( 20 ) 


(Demanding ipn( 7t/2) = ip* = 1 has fixed the constant 
to be C 2 = 1.) The radial magnetic field is plotted in 
Fig. ESI and shows even more pronounced bunching 
than in the monopolar case. 


where ipn( 7r /2) = ip* = 1 has fixed C 2 = 1. The radial 
component of the magnetic field at the horizon is plotted 
in Fig. 1 (c)[ which again shows the bunching. 


IV. SUMMARY 


C. Hyperbolic 

A third, “hyperbolic” solution was found in where 
the field is generated by a thin disk terminating at an in¬ 
ner radius b. The current and angular velocity functions 
satisfy our assumptions. For simplicity we will work in 
the limit b^> M, where they become]^] 


I(ip) = 4-7T Cl(lp) Ip, 

(17) 

n(ip)-n H . 

l-F VI ~ip 

(18) 


These expressions are a good approximation even when b 
corresponds to the innermost stable circular orbit. Per¬ 
forming the bunching calculation, we find 

I (w) 


We have shown that the high-spin bunching of field 
lines observed by TNM generalizes to arbitrary magnetic 
geometry under the assumption that the functional 
forms of the current I (ip) and field line angular velocity 
Cl(tp) both scale linearly with fin- Key to enabling 
this analytic argument was the existence of the Znajek 
condition at the horizon. This technique allows us 
to bypass the issue of having to solve the complete 
non-linear force-free problem to derive properties of 
the fields at the horizon. We gave a general argument 
and then illustrated the bunching with three specific 
magnetic geometries, which are plotted in Fig. [l] 
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(b) parabolic 




(c) hyperbolic 


FIG. 1. Bunching of field lines for radial, parabolic, and hyperbolic geometries. Diagrams of the field line geometry appear 
in the upper right. On the left, the dimensionless current T (black line) and angular velocity O (gray line) as a function of 
magnetic flux ip/ip* on the horizon ip < ip*, which we assume to hold at all spin. On the right, the radial magnetic field 
B r = deip/y/—g as a function of angle 9 < 7r/2 on the horizon for different values of the spin a = 0.1, 0.5, 0.9, 0.99,0.999, 0.9999 
(upwards on the left). 





















